Consensus Networks over Finite Fields 



F. Pasqualetti a , D. Borra b , F. Bullo a 

a Center for Control, Dynamical Systems, and Computation, University of California, Santa Barbara, USA 
b Dipartimento di Scienze Matematiche, Politecnico di Torino, Italy 

m 

y—i 

Abstract 

<N 

(— ( This work studies consensus strategies for networks of agents with limited memory, computation, and communication capa- 
bilities. We assume that agents can process only values from a finite alphabet, and we adopt the framework of finite fields, 
h— •} where the alphabet consists of the integers {0, . . . ,p — 1}, for some prime number p, and operations are performed modulo 
p. Thus, we define a new class of consensus dynamics, which can be exploited in certain applications such as pose estimation 
in capacity and memory constrained sensor networks. For consensus networks over finite fields, we provide necessary and 
sufficient conditions on the network topology and weights to ensure convergence. We show that consensus networks over finite 
i fields converge in finite time, a feature that can be hardly achieved over the field of real numbers. For the design of finite-field 
j^j consensus networks, we propose a general design method, with high computational complexity, and a network composition 
j^V rule to generate large consensus networks from smaller components. Finally, we discuss the application of finite-field consensus 
networks to distributed averaging and pose estimation in sensor networks. 
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^» 1 Introduction 

o 

Sensor and actuator networks have recently attracted interest from different research communities, and, in the 
last years, classic computation, control, and estimation problems have been reformulated to conform the distributed 
nature of these networked systems [1,2,3]. An important example is the so-called consensus problem, where members 
of a network aim to agree upon a parameter of interest via distributed computation [4] . Consensus algorithms have 
found broad applicability in many domains, including robotics [5], estimation [6], and parallel computation [7]. 



. . In this work we focus on the consensus problem for networks of agents with limited memory, computation, and 
^* communication capabilities. We assume that agents are capable of storing, processing, and transmitting exclusively 
k> elements from a finite and pre-specified alphabet, and we conveniently model this situation with the formalism of 
finite fields, where the alphabet consists of a set of integers, and operations are performed according to modular 
arithmetic [8]. We study linear consensus networks over finite fields where, at each time instant, each agent updates 
its state as a weighted combination over a finite field of its own value and those received from its neighbors. Besides 
consensus in capacity and memory constrained networks, our finite-field consensus method has broad applicability 
to problems in cooperative control, networked systems, and network coding, such as averaging, load balancing, 
and pose estimation from relative measurements. Additionally, the use of a finite alphabet for computation and 
communication makes our consensus method easily implementable and resilient to communication noise. 



Related work Consensus algorithms have been proposed for different network model, agents dynamics, and commu- 
nication schemes. Starting from the basic setup of time-independent network structure, broadcast and synchronous 
communication, and unlimited communication bandwidth, consensus algorithms have been proposed to cope, for 
instance, with time- varying topologies [9], gossip and asynchronous communication [10], and communication errors 
and link failures [11]. It has been shown that, under mild connectivity assumptions on the interaction graphs, a simple 
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linear iteration ensures consensus, where at each time agents update their state as a weighted average of the available 
neighboring states [12]. While most of these approaches assume the possibility of processing and transmitting real 
values, we consider the more realistic case of finite communication bandwidth, possibly due to digital communication 
channels, and memory constraints. Moreover, we will show that topological conditions ensuring consensus over real 
values and with real-valued communications are not sufficient for consensus over a finite field. 

A relevant body of literature deals with consensus over quantized communication channels, where values exchanged by 
the agents are quantized according to a predefined quantization scheme, and the proposed algorithms are resilient to 
quantization errors [13,14,15,16,17]. In these works, although the data exchanged among robots is quantized, agents 
perform computations over the field of real numbers. Thus, the consensus problem with quantized information is 
related to our problem, yet fundamentally different, since we allow agents to operate on a finite field. 

Logical consensus has been studied in [18] for the purpose of intruder or event detection. In logical consensus, a team 
of agents aims to coordinate their decisions via distributed computation as a function of a set of logical (boolean) 
events. By leveraging tools from cellular automata and convergence theory of finite-state iteration maps, the focus 
of [18] is on the design of a synthesis technique for logical consensus systems. The main differences between [18] and 
this paper are as follows. First, in logical consensus the agents state is a binary variable, while in our work it takes 
value in an arbitrary finite set. Second, in logical consensus agents are allowed to perform any logical operation, such 
as {and, or, not}, while our consensus algorithm makes use of only two mathematical operations, namely modular 
addition and multiplication. Third and finally, in logical consensus agents aim to agree upon a logical expression or 
compact sets [19], while finite- field consensus algorithms may be used to compute a (non-boolean) function of the 
initial states, such as the average. 

A distributed consensus algorithm with integer communication and computation is proposed and analyzed in [20]. 
With respect to this work, and similarly to [21], we make use of modular arithmetic, instead of standard arithmetic 
and, therefore, we define a novel and complementary class of consensus networks. The use of modular arithmetic 
is advantageous in several applications, such as pose estimation from relative measurements (Section 6). Finally, 
networks based on modular arithmetic are studied in [22] , in the context of system controllability and observability, 
in [23]. in the context of (linear) network coding, and in [24], in the more general context of finite dynamical systems. 

Contributions The contributions of this paper are fourfold. 

First, we propose the use of modular arithmetic to design consensus algorithms for networks of cooperative agents. 
Consensus networks over finite fields are distributed, require limited, in fact finite, memory, computation, and com- 
munication resources, and, as we show, they exhibit finite time convergence. Thus, finite-field consensus algorithms 
are suitable for capacity and memory constrained networks, and for applications with time constraints. 

Second, we thoroughly characterize convergence of consensus networks over finite fields. We provide necessary and 
sufficient constructive conditions on the network topology and weights to achieve consensus. For instance, we show 
that a network achieves consensus over a finite field if and only if the network matrix is row-stochastic over the finite 
field, and its characteristic polynomial is s rl_1 (s— 1). Additionally, we prove that the convergence time of finite-field 
consensus networks is bounded by the network cardinality, and that graph properties alone arc not sufficient to 
ensure finite-field consensus. Our analysis complements the classic literature on real-valued consensus networks. 

Third, we propose systematic methods to design consensus networks over finite fields. In particular, we derive a 
general design method, with high computational complexity, independent of the agents interaction graph, and a 
network composition rule based on graph products to generate large consensus networks from smaller components. 
We show that networks generated by our composition rule exhibit a specific structure, and maintain the convergence 
properties, including the convergence time, of the underlying components. Moreover, by using our general network 
design method, we compute finite-field consensus weights for some specific classes of interaction graphs, and we 
provide a lower bound on the number of networks achieving consensus as a function of the agents interaction graph 
and the field characteristic. 

Fourth and finally, we consider two applications in sensor networks, namely averaging and pose estimation from 
relative measurements. In the averaging problem agents aim to determine the average (over the field of real numbers) 
of their initial values. We show that, under a reasonable set of assumptions, the averaging problem can be solved 
distributively and in finite time by using a finite-field average consensus algorithm. In the pose estimation problem 
agents aim to estimate their orientation based on (local) relative measurements. For this problem we derive a 
distributed pose estimation algorithm based on finite-field average consensus, and we characterize its performance. 
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Paper organization In Section 2 we recall some preliminary notions on fields, linear algebra, and graph theory. 
Section 3 contains our setup and some preliminary results on consensus over finite fields. In Section 4 we provide 
necessary and sufficient conditions for consensus over finite fields, and we propose several illustrative examples. 
In Section 5 we describe methods to design consensus networks over finite fields. Finally, Section 6 contains two 
application scenarios, and Section 7 concludes the paper. 

2 Notation and Preliminary Concepts 

In this section we recall some definitions and properties of algebraic fields, linear algebra, and graph theory. We refer 
the interested reader to [8,25,26] for a comprehensive treatment of these subjects. 

A field F is a set of elements together with addition and multiplication operations, such that the following axioms 
hold: 

(Al) Closure under addition and multiplication, that is, for all a, b G F, both a + b G F and a ■ b G F; 
(A2) Associativity of addition and multiplication, that is, for all a,b,c G F, it holds a + (b + c) = (a + b) + c and 
a ■ (b ■ c) — (a ■ b) ■ c; 

(A3) Commutativity of addition and multiplication, that is, for all a, b G F, it holds a + b = b + a and a ■ b = b ■ a; 
(A4) Existence of additive and multiplicative identity elements, that is, for all a G F, there exist b, c € F such that 
a + b = a and a ■ c — a; 

(A5) Existence of additive and multiplicative inverse elements, that is, for all a G F, there exist b, c G F such that 

a + b = and a ■ c — 1, with a^O; 
(A6) Distributivity of multiplication over addition, that is, for all a, b, c G F, it holds a ■ (b + c) = (a ■ b) + (a • c). 

A field is finite if it contains a finite number of elements. A basic class of finite fields are the fields F p with characteristic 
p a prime number, which consists of the set of integers {0, . . . ,p — 1}, with addition and multiplication defined as 
in modular arithmetic, that is, by performing the operation in the set of integers Z, dividing by p, and taking the 
remainder. Unless specified differently, all the operations listed in this paper are performed in the field F p . 

Let a : F™ — >■ F™ be a linear map between the vector spaces of dimensions m and n, respectively, over the field F p . 
As a classical result in linear algebra, the map a can be represented by a matrix A with n rows and m columns, and 
elements from the field F p . The image and kernel of A are defined as 

Im(A) := {y G F™ : y = Ax, x G F™}, Ker(A) := {x G F™ : Ax = 0}, 

where additions and multiplications are performed modulo p. Analogously, the pre-image of a set of vectors V C F p 
through A is the set 

A _1 (V) := {x G F™ : v = Ax, for all v G V}. 
For a matrix A G F p x ", let cf v {A) denote the set of eigenvalues of A in the field F p , that is, 

o- p (A) = {A G F p : det(AJ - A) = 0}. 

Analogously, let F p [s] denote the set of polynomials with coefficients in F p , and let Pa G F p [s] denote the char- 
acteristic polynomial of A over F p PH Notice that the cardinality |cr p (A)| may be strictly smaller than the matrix 
dimension n; in other words, finite fields are not algebraically closed. 

We conclude this section with some standard graph definitions. A directed graph Q = (V,£) consists of a set of 
vertices V and a set of edges £ C V x V. An edge (v,w) G £ is directed from vertex w to vertex v. For a vertex 
v G V, the set of in-neighbors of v is defined as A^ n = {w G V : (v,w) G £}, and the set of out-neighbors as 
Af° ut = {w G V : (w, v) G £}. The in-degree of v G V equals |A/^ n |, whereas the out-degree of v G V equals |A/"° ut |. 
A path in G is a subgraph P — ({ui, . . . , ffc+i}, {ei, . . . , e/c}) such that 7^ Vj for all i 7^ j, and = (u^, for 



1 Since the characteristic polynomial Pa(s) G R[s] contains only integer coefficients for any A G Fp Xn , the characteristic 
polynomial Pa(s) G F p [s] is J]™ =0 mod(ci,p)s\ where mod(-) is the modulus function, and Ci is the i-th coefficient of Pa(s). 
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each i G {1, . . . , k}. We say that the path starts at vi and ends at Vk+i, and, at times, we identify a path by its 
vertex sequence v±, . . . , Vk+i- A cycle or closed path is a path in which the first and last vertex in the sequence are 
the same, i.e., V\ = Vk+i- The length of a path (resp. cycle) equals the number of edges in the path (resp. cycle). 
A directed graph is strongly (resp. weakly) connected if there exists a directed (resp. undirected) path between any 
two vertices. Two subgraphs of the same graph are disjoint if they have no common vertices. A root (resp. globally 
reachable node) is a vertex v from which (resp. to which) there exists a directed path to (resp. from) every vertex in 
the graph, including v itself. Finally, a directed graph is aperiodic if there is no integer greater than one that divides 
the length of every cycle of the graph. 

3 Models of Consensus Networks over Finite Fields 

Consider a set of n G N>o agents and a finite field F p , for some prime number p. Let the agents interaction be 
described by the directed graph Q = (V, £), where i E V denotes the i-th agent, with V = {1, . . . , n}, and G £ 
if there is a directed edge from agent j to agent i (agent i senses agent j, or, equivalently, the behavior of agent 
j affects agent i). We assume that each agent is able to manipulate and transmit values from the finite field F p 
according to a pre-specified protocol. We focus on distributed protocols in which (i) each agent i is associated with 
a state Xi G F p , and (ii) each agent updates its state as a weighted combination of the states of its in- neighbors A^ 1 ". 
Let ciij G F p be the weight associated with the edge and let A = [(%■], A G F p xn , be the weighted adjacency 

matrix of Q, or simply network matrix, where = whenever G" £■ Let x : N>o — > F p be the vector containing 
the agents states. Then the evolution of the network state x over time is described by the iteration (or network) 

x{t + l) = Ax{t), (1) 

where all operations are performed in the field F p . 

The transition graph associated with the iteration (1) over F p is defined as Qa — {Va, £a), where, Va — {v : v G F p } 
and, for v i} Vj G Va, the edge (vi,Vj) G £a if and only if Vj = Av; t . It should be observed that the transition graph 
contains p n vertices, and that each vertex has unitary out-degree. Moreover, it can be shown that the transition graph 
is composed of disjoint weakly-connected subgraphs, and that each subgraph contains exactly one cycle, possibly of 
unitary length [27]. Finally, each disjoint subgraph contains a globally reachable node. This particular structure of 
the transition graph will be used to derive certain results on finite-field consensus. Examples of transition graphs 
are given below in Fig. 1 and Fig. 2. 

We say that the iteration (1) (or simply the network matrix A) over a finite field achieves 

(i) asymptotic consensus, if for all initial states x(Q) G F p it holds lim^^ x(t) = at, with a G F p and 
1 = [1...1] T ; 

(ii) finite-time consensus, if for all initial states x(0) G F p there exists a finite time T G N such that x(T) = 
x(T + t) = at for all r G N, with a G F p and 1 = [1 . . . 1] T . 

Consensus networks with real- valued weights and states have been extensively studied [12,28,4]. In this work we show 
that real-valued consensus networks and finite-field consensus networks exhibit different features, and particular care 
needs to be taken to ensure the desired properties over finite fields. It is clear from the above definitions that finite- 
time consensus implies asymptotic consensus. We next show that the converse is also true. 

Theorem 3.1 (Asymptotic consensus implies finite-time consensus) The iteration (1) over the field F p 
achieves asymptotic consensus only if it achieves finite-time consensus. 



PROOF: Let Qa = (Va,£a) be the transition graph associated with the iteration (1). Notice that the state 
trajectory x of (1) coincides with a path on Qa starting from the vertex vo = x(0). Let C C Va be the set of 
consensus vertices, that is, C = {v : v G Va,v = at, a G F p }. Suppose that the iteration (1) achieves consensus 
on the value v c G C. Since the vertex set Va is finite, the distance between vq and v c is also finite. Consequently, a 
consensus vertex is reached with a path on Q of finite length, that is, a finite number of iterations in (1) are sufficient 
to achieve consensus. □ 

Following Theorem 3.1, iterations over finite fields either achieve consensus in finite time, or they are not convergent. 
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In view of this result, in what follows we simply use consensus instead of finite-time and asymptotic consensus. 
Differently from finite-field consensus networks, consensus networks over the field of real numbers usually converge 
asymptotically. An exception is constituted by the class of de Bruijn graphs, which have been shown to yield finite- 
time consensus over the field of real numbers [29]. On the other hand, de Bruijn graphs rely on a specific interaction 
graph, while finite-field consensus networks include a much broader class of interaction graphs. We conclude this 
section with a simple result. A matrix A over the field F p is nil-potent if A n = and is row- stochastic if At = 1. 

Lemma 3.2 (Finite- field consensus matrices) Consider the iteration (1) over the field F p . If consensus is 
achieved, then A is either nil-potent or row-stochastic. 



PROOF: Since A achieves consensus, it follows from Theorem 3.1 that A t x{Q) = A t+1 x{0) = at for some aeF n , 
for all x(0), and for all t > T, T e N. Then Aat = at, from which we conclude that either At = 1 (A is row- 
stochastic) or a — for all initial states a;(0) (A is nilpotent). □ 

As for the case of real- valued consensus, we limit our attention to row-stochastic network matrices. Although con- 
sensus is trivially achieved whenever the network matrix is nilpotent, this case is of limited interest because the 
consensus value is the origin independently of the agents initial states. 

4 Analysis of Consensus Networks over Finite Fields 

Conditions for consensus in real- valued networks have been deeply investigated in the last years [12,28,4]. For instance, 
sufficient conditions ensuring real-valued consensus are that the network matrix A is row-stochastic and that the 
associated directed graph is strongly connected and aperiodic. The following example shows that graph-theoretic 
properties are not sufficient for an iteration over a finite field to achieve consensus. 

Example 1 (Graph properties are not sufficient for finite-field consensus) Consider a fully connected 
network with three agents over the field F 3 . Consider the network matrices 
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Notice that A\, A 2 , and A 3 are row- stochastic and their interaction graph is fully connected. It can be verified that 
over the field F 3 only the network matrix A\ achieves consensus, while A 2 and A 3 exhibit oscillatory dynamics for 
certain initial conditions. An example of oscillatory dynamics generated by A 3 is reported in Table 1. 

Table 1 

Sample state trajectory for the matrix A A in Example 1. 
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□ 

As shown in Example 1, graph properties of the network matrix are not sufficient to guarantee consensus for iterations 
over a finite field (in short, finite-field consensus). Indeed, although the considered network matrices feature the 
same connectivity properties, only one of them achieves finite-field consensus. In what follows we provide finite-field 
consensus conditions based on the algebraic properties of the network matrix, and on the topological properties of 
its transition graph. 

The dynamic behavior of an iteration over a finite field is entirely described by its associated transition graph. The 
next theorem provides a necessary and sufficient condition for finite-field consensus based on the transition graph. 

Theorem 4.1 (Transition graph of a consensus network) Consider the iteration (1) over the field F p with 
row-stochastic matrix A, and let Ga — (Va,£a) be its associated transition graph. The following statements are 
equivalent: 
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Fig. 1. Transition graph Qa x associated with the matrix A\ G Fg X3 in Example 1. Since Qa x contains exactly 3 cycles 
corresponding to the self-loops around the consensus vertices, the network matrix A\ achieves consensus (see Theorem 4.1). 




Fig. 2. Transition graph Qa 3 associated with the matrix A$ G Fg X3 in Example 1. Since Qa 3 contains more than 3 cycles, the 
network matrix A3 does not achieve consensus (see Theorem 4.1). The oscillatory state trajectory in Table 1 corresponds to 
the bottom right cycle in this figure. 

(i) the iteration (1) achieves consensus, and 

(ii) the transition graph Qa contains exactly p cycles, corresponding to the unitary cycles around the vertices at, 
a g {0,...,p- 1}. 

Example 2 (Transition graph of a consensus network) The transition graphs associated with the matrices 
A\ and A3 in Example 1 over the field F3 are reported in Fig. 1 and Fig. 2, respectively. As previously discussed, 
and as predicted by Theorem 4-1, the matrix A\ achieves consensus, while the matrix A3 does not. □ 



PROOF of Theorem 4.1: 

(i) => (ii) Since the iteration achieves consensus, it follows from Lemma 3.2 that At = t. Hence, the transition 
graph contains p unitary cycles corresponding to the vertices at <G Va, with a e F p . Suppose by contradiction that 
there exists an additional cycle C, and notice that the vertices at, with a € F p , cannot be contained in C since 
the out-degree of each vertex in the transition graph is exactly one (the transition graph is determined by the linear 
map A). Thus, there exists a trajectory along C that does not converge to consensus, which contradicts the initial 
hypothesis. 

(ii) =>■ (i) Notice that a state trajectory of the iteration (1) is in bijective correspondence with a path on the tran- 
sition graph Q A . Suppose that transition graph Q A contains exactly p unitary cycles located at the vertices at e Va, 
with a g F p . Then, since each vertex in the transition graph has unitary out-degree, every (sufficiently long) path in 
Qa eventually reaches one of the cycles, and, consequently, every state trajectory converges to a consensus state. □ 



Theorem 4.1 provides a necessary and sufficient condition for finite-field consensus based on the transition graph. 
From condition (ii) in Theorem 4.1 and the fact that each vertex in the transition graph has unitary out-degree, 
we also note that the transition graph of a consensus matrix is composed of p disjoint weakly-connected subgraphs. 
Moreover, by means of [27, Proposition 3.4], it can be shown that disjoint subgraphs have the same graph topology. 

A verification of the convergence condition in Theorem 4.1 may be prohibitive for large networks, because the size of 
the transition graph grows exponentially with the number of agents in the network (the transition graph contains p n 
vertices and p n edges, since each vertex has unitary out-degree). In what follows we shall derive consensus conditions 
based on the network matrix instead of its transition graph. Consider the inverse recursion 

Si+^A-^Sl), (2) 
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where S^ C F™ for all times t and <S° = {al}. Notice that the inverse recursion defines a sequence of sets, and that 
the set <S^ contains the initial states converging to the consensus value at in at most t iterations. We say that the 
recursion (2) is convergent with limiting set S a if there exists T < n satisfying S a — S T = <?J +1 . The following 
theorem exploits the inverse recursion and the structure of the transition graph to characterize finite-field consensus. 

Theorem 4.2 (Recursion subspaces of a consensus network) Consider the iteration (1) over the field F p 
with row-stochastic matrix A. The following statements are equivalent: 

(i) the iteration (1) achieves consensus, 

(ii) there exists a£F p such that the recursion (2) is convergent and the limiting set S a satisfies \S a \ — P n ~ l , and 
(Hi) for all a e F p the recursion (2) is convergent and each limiting set S a satisfies \S a \ = p n ~ x . 

Example 3 (Inverse recursion for a finite-field consensus network) For the matrix A\ e F 3x3 in Example 
1, the set Si generated by the inverse recursion (2) is 
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Because |<Si| = 3 2 , the network matrix A± achieves consensus due to Theorem 4-2. Instead, for the network matrix 
A 3 e Fg X3 in Example 1, the inverse recursion yields Si — {t}, so that A 3 does not achieve consensus. □ 



PROOF of Theorem 4.2: Consider the transition graph Qa — (V, £), and define the reverse graph Ga = (V, £), 
where G £ if and only if (j, i) G £. Notice that the recursion (2) is convergent if and only if Ga contains no 
cycle of length greater than 1 reachable from at. Recall from [27, Theorem 1] that Ga (resp. Ga) is obtained as the 
graph product of a tree by a set of cycles. Hence, the graph Ga (resp. Ga) is composed of disjoint weakly-conncctcd 
subgraphs, and disjoint subgraphs have the same structure. From this argument we conclude that (ii) and (iii) are 
equivalent. 

(i) (ii) Since A achieves consensus, the graph Ga contains exactly p unitary cycles corresponding to the 
consensus vertices (see Theorem 4.1 and Fig. 1). By [27, Theorem 1], the above reasoning, and the fact that A 
achieves consensus, it follows that |<So| H + = p n , and that |<S Q | = p n ~ x for all a e F p . 

(ii) =>■ (i) Since At = 1, the transition graph contains p cycles of unitary length located at the consensus vertices. 
Let the recursion (2) be convergent for some a € F p . From [27, Theorem 1], the graph Ga contains p identical, 
disjoint, weakly-connected subgraphs, each one terminating in a consensus vertex. Since \S a \ = p n ~ x , it follows that 
consensus is achieved from p n states (every initial state), which concludes the proof. □ 



According to Theorem 4.2, the convergence of the network (1) can be determined by iterating the inverse recursion 
(2) for some a e F p . This computation does not require analyzing the transition graph. Our last and most explicit 
condition for finite-field consensus is based upon the characteristic polynomial of the network matrix (computed 
over the finite field). 

Theorem 4.3 ( Characteristic polynomial of a consensus network) Consider the iteration (1) over the finite 
field F p with row- stochastic matrix A. The following statements are equivalent: 

(i) the iteration (1) achieves consensus, and 
(ii) P A (s) = s n - 1 (s-l). 

Example 4 (Characteristic polynomial of consensus matrices) Consider the network matrices in Example 
1 over the field F 3 . It can be verified that 

PaM = s 2 (s - 1), Pa 2 (s) - s(s 2 - 2s + 1), and P A3 (s) = s 3 - 1. 

As predicted by our previous analysis and by Theorem 4-3, only the network matrix Ai achieves consensus. □ 
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The proof of this theorem is postponsed to the Appendix. Theorem 4.3 is cquivalcntly restated as follows: A achieves 
finite-field consensus if and only if a p (A) = {1,0, ...,0}. In other words, the eigenvalues of a finite-field matrix 
achieving consensus are all contained in the considered finite field and, consequently, every finite-field matrix achiev- 
ing consensus can be represented in Jordan canonical form via a similarity transformation; see [30] and [31, Theorem 
3.5]. We conclude this section by characterizing the convergence value of a finite-field consensus network. 

Theorem 4.4 (Finite-field consensus time and value) Consider the iteration (1) over the finite field F p with 
row-stochastic matrix A and with initial state x(0) £ F™. Assume the iteration (1) achieves consensus. Let T < n 
denote the dimension of the largest Jordan block associated with the eigenvalue 0. Let tt £ F™ be the unique eigenvector 
satisfying irA = tt and nt = 1 . Then 

A T = Itt, 

so that consensus is achieved at the value irx(0) after T iterations. Moreover, the i-th component of ir is nonzero 
only if the i-th vertex of the directed graph associated with A is a root. 



PROOF: Since A achieves consensus, we have o- p (A) = {1,0,..., 0}, and A admits a Jordan canonical form 
J A = V~ 1 AV over F p . Moreover, the matrix A converges in T < n iterations. The next part of the proof follows the 
reasoning in [32, Theorem 3]. Let the first column of T be 1, and let the matrix J a have only zero elements, except 
for the unitary entry in position (1, 1). Since V^ 1 A = JaV^ 1 , the first row of V^ 1 , say w, satisfies irA — tt. Then 
A T = VJ\V~ X = Itt. Since At = 1, it follows that 1 = A T t = tirt, and consequently irt = 1. To show the last 
statement, let Q be the directed graph associated with A, and let i be a vertex of Q. Assume that i is not a root of Q, 
and let the initial state x(0) be all zeros, except for the i-th component. Since i is not a root, there exists a node j 
which is not reachable from i, and, consequently, the value of the j-th agent is not affected by the i-th agent. Since A 
achieves consensus for all initial states, the j-th entry of 1ttx(0) needs to be zero, from which the statement follows. □ 



Observe that Theorem 4.4 is not a direct consequence of the theory of non-negative matrices over the field of real 
numbers [30]. In fact, if regarded as a real- valued matrix, a finite-field consensus matrix is generally unstable. 

5 Design of Consensus Networks over Finite Fields 

In Section 4 we characterize the convergence of consensus networks over finite fields. With respect to consensus 
networks over the field of real numbers, finite-field consensus networks require less computational effort and commu- 
nication bandwidth, and they converge in a finite number of iterations. On the other hand, convergence conditions 
for finite-field consensus networks depend on the numerical entries of the network matrix (Theorem 4.3), and not 
only on the connectivity properties of the underlying graph as in the case of consensus networks with real values. 
For this reason, the design of finite-field consensus networks deserves particular attention. In this section we describe 
methods to design finite-field consensus networks, and we discuss their limitations. First, we elaborate on Theorem 
4.3 to describe a general design method, and we propose a solution for particular interaction graphs. Then we present 
a scalable composition rule to generate finite-field consensus networks from smaller consensus components. 

5. 1 Network design via characteristic polynomial 

The objective of this section is to design a finite-field consensus matrix A whose sparsity patter is compatible with 
a given agents interaction graph Q = (V, £), that is, to design the network matrix A = [ay], ay £ F p , where p is a 
given prime number, and ay ^ only if (i, j) is an edge of Q . Recall from Theorem 4.3 that the network matrix A 
achieves consensus if and only if its characteristic polynomial is Pa{s) = s n (s— 1) and At = t. Since two polynomials 
are equal if and only if all coefficients are equal, the entries of A can be determined by simultaneously solving the 
following equations: 

' A £ F" x ", with ay = if (i, j) <£ E, 
1 = A1, " 

1 =-c(A,n-l), ( ' 

(0 =c(A,j), je{0,...,n-2}, 
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where c(A,d) is the coefficient of the monomial of degree d in the (parametric) characteristic polynomial Pa- 



Example 5 (Design for networks with 2 agents) Let A 

(an +022)5 + 011022 — ai2a2i- It follows from (3) tha, 
0-12 + &21 = I, and, consequently, if and only if A = 



c F 2x2 
13 ^ p > 



an oia 

a 21 a 22 

(011+022)5 + 011022 — ai2Q-2i- It follows from (3) that A achieves consensus if and only if an — 1 — 012, 022 

1 — a a 



1 



a 



for some a G F 



and notice that Pa(s) = s 2 ~ 

1-021, 

□ 



Notice that the system of equations (3) contains non-linear, multivariate, polynomial equations, where the unknown 
variables are the entries of the network matrix A, and that a solution is required over the finite field F p . The 
problem of solving systems of multivariate polynomial equations over finite fields is NP-hard [33,34], and it is one 
of the important research problems in cryptography and information security [35]. Besides enumerating all possible 
solution candidates]^] several solution techniques have been proposed over the last years, see for instance [36,37]. 
We next provide a condition on the agents interaction graph for the existence of a finite-field network achieving 
consensus. 

Theorem 5.1 (Existence of finite-field consensus matrices) Let Q = (V,£) be the directed agents interaction 
graph, with \V\ — n and \£\ — m. Assume that Q contains a root, and that m > " 2 1 '" . Then the system of equations 

2m-„ 2 -n 

(3) admits N > p 2 solutions over the field ¥ p , where N is divisible by the characteristic p. In other words, 
there exist N network matrices A = [ay] achieving consensus, with atj € ]?™ xn an & a ^ = if (i,j) ^ £ . 



PROOF: Let c(A, d) be the coefficient of the monomial of degree d in the parametric characteristic polyno- 
mial Pa, with A = [a^] and ay = if (i,j) (I £■ Define the polynomial fi € F p [oi :) ] with (i,j) G £ as: 
ft = /i(fflij) = Xy=i a ij - 1 for i e {1, ... , n}, f n+ i = 1 + c(A, n - 1), and /, = c(A, 2n - i) for i e {n + 2, . . . , 2n}. 
Notice that a solution to (3), and hence a finite-field consensus matrix compatible with Q, can be computed by 
simultaneously solving the equations fi — for i G {1, . . . , 2n}. Observe that deg(/;) = 1 for i E {1, ... ,71}, and 

deg(/») < i - n for % € {n+ 1, . . . , 2n}. Then deg(/i) < n + < . Since m > by assumption, we 

conclude from [8, Theorem 6.8] that the number of simultaneous solutions N to the equations fi = 0, i € {1, . . . , 2n}, 

2 m — n 2 —n 

in the field F p is divisible by the characteristic p. To show that N > p 5 , we next construct a network matrix 
A achieving consensus and compatible with the interaction graph Q, and then employ [8, Theorem 6.11]. Let v be a 
root of G, and let S = (V,£s) be a rooted spanning tree of Q with root v, and with (v,v) £ £s and £ £s for 
i =/= v. Relabel the vertices V according to their distance from the root v. Define the matrix A = [ay], where ay = 1 
if (i,j) £ £s, and ay = otherwise. Notice that A € F™ xn is triangular and row-stochastic, and that its diagonal 
elements are {1, 0, . . . , 0}. It follows from Theorem 4.3 that A achieves consensus. □ 



In view of the existence Theorem (5.1), consensus over finite fields is possible on a fairly broad class of interaction 
graphs. A complete characterization of all the interaction topologies yielding consensus over finite fields is beyond 
the scope of this work, and it is left as the subject of future research. We conclude this section with a remark. 

Remark 1 (Network design for fully connected graphs) Let the agents interaction graph Q ~ (V,£) be fully 
connected, that is G £ for all i,j G {1, . . . ,n\. Let v £ F p x ™ be any vector satisfying vl = 1. Then the network 
matrix A = [v T ■ ■ ■ v J ] T achieves consensus over F p ; see for instance the matrix A\ in Example 1. To see this, let 

^orth ^ F p x ^" 1 ' be any full column rank matrix satisfying vt mt ^ = 0. Since At = 1, and j4l orth = 0, we conclude 
that CFp(A) = {1,0,..., 0}, with |cr p (A)| = n, and the claimed statement follows from Theorem 4-3. □ 



5.2 Network design via network composition 

In this section we use the concept of graph products to generate finite-field consensus networks from smaller consensus 
components. For two matrices A e F™ x " and B e FJ 1 *"", let A <£> B e F" mxnm denote their Kronecker product, 



2 If m is the number of free entries in A, a brute-force solution to (3) requires computing all p m possible matrices A compatible 
with the interaction graph Q. 
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o 








Fig. 3. Agents interaction graphs for the matrices A (left) and — Ac 
Notice that the interaction graphs of A and A^ are self-similar [38]. 



'A (right) in Example 6. Self-loops have been omitted. 



where [30] 



A® B = 



ai\B CI22B 
a n \B a n 2B 



a\ n B 

®n.T}.B 



Theorem 5.2 (Finite- field consensus via Kronecker product) Consider the network matrices A e F™ 
B e jprnxm^ an d assume that A and B achieve consensus. Then the network matrix A ® B G 
consensus. 



and 



¥ nmxnm achieves 



Example 6 (Finite- field consensus network via Kronecker product) Consider the network matrix 



A = 



9 3 
1 9 2 
7 5 



over the field Fn. It can be verified that A achieves consensus. By Theorem 5.2 the network matrix 



Ay = A ® A = 



4 5 5 
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achieves consensus over Fn. In fact, it can be verified that Pa^ — s 8 (s — 1) and A^l = 1. Moreover, the network 
matrices A and A k have the same convergence speed. The interaction graph of A and A k are reported in Fig. 3. □ 
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PROOF of Theorem 5.2: Wc first show that (A ® B) is row-stochastic. Let 1„ be the vector of all ones of 



dimension n. Since A and B are row-stochastic over F p , we have 



{A®B)t r 



auB a 12 B 
ai\B a,2iB 

a n \B a n 2B 



a\ n B 

CI n r i B 



ho 



Bt r 



En 
, = l a 
Ej=l a 2,jB± Ti 



Ej=l a l,j : "-m 



Le.; 



=1 ^n,j^-m. 



Following the same argument as in [30], let J a = P 1 AP and Jb = Q X PQ be the canonical Jordan form of A 
and B, respectively [31, Theorem 3.5]. Then J A ® Jb = (P _1 AP) <g> (Q^BQ) = {p- 1 ® Q'^A ® B)(P ® Q) = 
(P®Q)~ 1 (A®B)(P®Q), so that A®B and Ja®Jb are similar. Thus the eigenvalues of A®B are the same as those 
of J a® Jb, and, because J a and Jb are upper triangular with eigenvalues Aj and /i, on the diagonal, we conclude 
that Ja®Jb is also upper triangular with diagonal entries, and eigenvalues, Xifij. Finally, since a p (A) = {1, 0, . . . , 0} 
and <J P (B) = {1,0,..., 0}, the statement follows from Theorem 4.3. □ 



Following Theorem 5.2, finite-field consensus networks can be constructed by composing smaller components. We 
refer the interested reader to [39,38] for a comprehensive discussion of graphs generated via Kronecker product of 
adjacency matrices. Regarding the convergence speed of finite-field consensus networks generated via Kronecker 
products, the following holds. Let A^ = A\ ® ■ ■ ■ ® A m , and let Sj be the number of iterations needed for convergence 
of the network Aj, i e {1,...,™}. Due to [40, Theorem 4.3.17], the consensus network A^ converges exactly in 
maxjs!, . . . , s m } iterations. In other words, the consensus network A^ is as fast as the slowest of its components. 



6 Application to Average Consensus and Distributed Pose Estimation 



This section contains two application scenarios for finite-field consensus networks. In Section 6.1 we develop a finite- 
time averaging algorithm based on finite-field consensus networks. Instead, in Section 6 we use finite-field networks 
to distributively estimate the orientation of a sensor network given relative measurements. 



6.1 Finite-time average consensus 



Given a sensor network, let xq € F™ be the vector containing the agents initial states. Let xr = I^xq/u etbe the 
average of the agents initial states over the field of real numbers. The average of the agents initial states over the 
field F p can be computed by means of Fermat's little theorem [41] as x^ = n p ~ 2 t T XQ e Fp, where we assume n ^ kp, 
with k e N, for the inverse of n over F p to exist. In what follows, first we show how to compute the average Xf 
by means of finite-field consensus networks. Then we describe conditions that allow to recover the average xr from 
the knowledge of x^ and the total number of agents. We remark that distributed algorithms for average consensus 
are useful in several applications, including distributed parameter estimation [6], distributed hypothesis testing [42], 
robotic coordination [2], and for the understanding of opinion dynamics [43]. 

We say that the iteration (1) over the field F p achieves average consensus if it achieves consensus, and the consensus 
value is n p ~ 2 t T xo for every initial state xq. A condition for finite-field average consensus is given in the next theorem. 

Theorem 6.1 (Finite-field average consensus) Consider the iteration (1) over the field F p with row- stochastic 
matrix A. Assume that the field characteristic satisfies n ^ kp for all k € N. The following statements are equivalent: 



(i) the iteration (1) achieves average consensus, and 
(it) Pa(s) = s" _1 (s - 1), and t T A '= 1 T . 

Example 7 (An example of finite-field average consensus) Consider the network matrix 



A = 



2 3 1 
2 4 
2 4 
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Fig. 4. A subgraph of the transition graph associated with the network matrix A in Example 7. Notice that the sum of the 
initial states is maintained, and thus average consensus is achieved. 

over the field F5. It can be verified that At = t, t T A = 1 T , and Pa = s 2 (s — 1). By Theorem 6.1 the network matrix 
A achieves average consensus over F 5 . In Fig. 4 we show a subgraph of the transition graph associated with A. □ 



PROOF of Theorem 6.1: 

(i) ==>■ (ii) Since the iteration achieves consensus, it follows from Theorem 4.4 that A n = tir, where ir satisfies 
ttA = 7r. Because A achieves average consensus, it needs to be tir = n p ~ 2 tt T . Then tt = n p ~ 2 t T , and t T A = t T . 

(ii) (i) Since Pa{s) = s n ^ 1 (s — 1) and At = 1, it follows from Theorem 4.3 that the network achieves 
consensus. Notice that t A = 1 T implies that t T x(t) — l T x(0) at all times time t. Let a be the consensus value, 
and notice that na = t T x(0). To conclude the proof, a = n p ~ 2 t ' x(0), and the network achieves average consensus. □ 



Theorem 6.1 provides a necessary and sufficient condition for a network with n ^ kp agents to achieve average 
consensus over F p . The condition n 7^ kp is actually necessary for average consensus. In other words, if n = kp for 
some k e N, then there exists no network matrix satisfying all conditions in Theorem 6.1 and, therefore, average 
consensus cannot be achieved. To see this, let x(Q) be the network initial state, with t T x(0) 7^ 0, and assume by 
contradiction that a is the corresponding consensus value. Since t T x(t) = l T a;(0) at all times t, it needs to be 
na = t T x(0). Then = na = l T x(0) 7^ 0, since kp and are in fact the same element in F p . 

Suppose now that the average xw has been computed, and that each agent knows the total number of agents, the 
field characteristic, and its own initial state. With these assumptions, it is generally not possible to recover the 
average xr. To see this, consider the case n = 3, p = 5, and the initial conditions X\ = [2 2 2] T and X2 — [0 1] T . 
Over the field of real numbers we have £i ; r = t T x\/n — 2 and X2,r = t J X2/n = 1/3. Over the field F p , instead, 
ill = n p ~ 2 t T xi = 2 and o^f = n p ~ 2 t T X2 = 2. Since Xij = x 2 ,w and a; m 7^ X2.R, it is not possible to recover the 
average value over the field of real numbers from the average over a finite field and knowledge of network cardinality 
and parameters. In the next theorem we present a sufficient condition to recover the desired real-valued average. 

Theorem 6.2 (Average computation) Let x e F™ 7 let x R = t T x a /n e R, and let x ¥ = n p ~ 2 t T x . If the field 
characteristic satisfies n||xo||co < P, then x^ = mod(nxr,p)/n. 

PROOF: The statement follows from the relation 

mod(n Xf , p) = mod(ti p_1 1 1 x a , p) = t T x , 
where the last equality holds because mod(n p_1 ,p) = 1, and n||x ||oo < P- □ 



Our finite-time average consensus algorithm follows from Theorems 6.1 and 6.2, and it is reported in Algorithm 1. 

We conclude this part by noticing that the condition 1 1 ccq 1 1 oo ^ p in Theorem 6.2 is not restrictive. In fact, the field 
characteristic p is a design parameter, and, in general, it can be chosen to satisfy the above condition as long as the 
network cardinality n is known and a bound on the agents initial state is known. 
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Algorithm 1 Distributed average computation ( agent i ) 

Input: Field characteristic p, Initial state x»(0) G F p , Neighbors set A^ m and JV° ut , Weights a tJ G F p for all j G Af™, 
Number of iterations T (set T = n otherwise), Number of agents n; 

Require: p is a prime number, A = [oij] achieves average consensus over F p , nmax{ii(0), . . . ,x n (0)} < p; 
Output: Average of initial states a; avo = Y^j=i x A^)/ n £ 

for i = 0, . . . , T do 

Transmit x l (t) to A/? ut ; 
Receive aij(t) from A/J n ; 
Set Xi(t+ 1) = I] jeA r|" aijXj{t); 
end for 

return a; a ve = mod(n Xi(T) , p) / n; 



6.2 Pose estimation from relative measurements 



In this section we use our previous analysis to calibrate the orientation of a network of cameras. This problem has 
been previously considered in [44,45,46] as a distributed estimation problem over 50(2). With respect to the existing 
literature, we let the measurements and the orientations take value in a pre-specified finite field, and we develop an 
estimation algorithm with performance guarantees based on modular arithmetic. 

A camera network is modeled by an undirected graph Q = (V, £), where each vertex is associated with a camera. 
Let n = \V\ and m — \£\. Let 0, : N>o O p be the orientation of the z-th camera as a function of time, where, for 
some prime number p, 




We refer to p as to discretization accuracy. For notational convenience, we define the directed graph = (Vd,£d) 
associated with the camera network Q, where Vd = V, and G £d if and only if G £ and i < j. For each 
(i,j) G £d, let r/ij G O p be the relative measurement between camera i and camera j, that is r/ij = 9i — 8j. Let 9 be 
the vector of the cameras orientations, and let rj be the vector of relative measurements. Assign an arbitrary ordering 
to the edges £d, and define the incidence matrix B G F™ x ™ of £d by specifying the fc-row of B corresponding to the 
edge as 

ri, ]f£ = i, 

b kl =\-l, ii£ = j, (5) 
[fj, otherwise. 

Observe that B9 = v. We consider the following problem. 

Problem 1 (Pose estimation from relative measurements) Let Q be a camera network, and let be its 

associated directed graph. Let p be the discretization accuracy, let B be the incidence matrix ofQ^, and let n G 0™ 
be the vector of relative measurements. Determine a set of cameras orientations 9 G 0™ satisfying B9 = r\ over F p . 



Notice that, if the incidence matrix B and the relative measurements n are available to some camera or central unit, 
then Problem 1 requires the solution of system of linear equations. Instead, we propose an algorithm that requires 
each camera to have access to local relative measurements and to communicate with its immediate neighbors. Our 
distributed pose estimation algorithm is in Algorithm 2. 

Because cameras transmit and operate only on values in the finite field F p , we argue that Algorithm 2 is suitable 
for agents with limited capabilities, and it is robust to transmission noise. In the next theorem we analyze the 
convergence of Algorithm 2. 

Theorem 6.3 (Convergence of Algorithm 2 with perfect measurements) Let Q = (V, £) be a camera net- 
work, and let be its associated directed graph. Let n G O™ be the vector of relative measurements, and let B be 
the incidence matrix ofQ d . If n G lm(B) and A achieves average consensus over ¥ p , then 
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Algorithm 2 Distributed pose estimation ( camera i ) 



Input: Discretization accuracy p, Initial pose 9i(0) € O p , Neighbors set Mr and M° ut , Weights aij G F" xn for all j G Mr, 
Number of iterations T (set T = n otherwise), Relative measurements r)ij G O v for all j G A/j ln ; 
Require: p is a prime number, A = [ojj] achieves average consensus over F p ; 
Output: Orientation 9i compatible with measurements n — [rjij]; 

for t = 0, . . . , T do 

Transmit Xi (t) = to M° ut ; 

Receive Xj(t) = from A/" 4 m ; 

Update orientation 0i(t) as: 

0;(t + l) = a: i (t + l)y. (6b) 

end for 

return Orientation 0<; 



(%) Algorithm 2 converges in finite time, that is, 9 := 9{T) = 9{T + r) for some 9 g (9™ , T < n, and for all r g N, 
and 

(mJ £/ie /ina/ network orientation is compatible with the relative measurements, that is, B9 = n. 



PROOF: Consider the update law (6a), and notice that it can be written as x(t + 1) = Ax(t) + LBv, where 
M = q v f or some vector v g F™ (y g Im(_B) by assumption), L g F™*" 1 , and the fc-th column of i corresponding 
to the edge € £d is specified as 



I' Ik 



Observe that 



(7) 



otherwise. 



{LBh 



-ay, ifjeA^ in , 
0, otherwise. 



so that LBt = (asymmetric Laplacian matrix of Q [2]). Notice that x(t) = A t x(0) + J2 t T J a A T LBv. Since A achieves 
average consensus, we have A T = n p ~ 2 tt T for some T < n. Thus, A f LB = n p ~ 2 tt T LB for all t > T. We now show 
that 1 T LB = 0, from which statement (i) follows. Since A achieves average consensus, it follows from Theorem 6.1 
that At = 1 and t T A = t T . Hence, for each node k g V, Y^j=i a kj — Ej=i a jk, and, consequently, t T LB = 0. 



Let i be a fixed point of (6a), that is, (/ — A)x = LBv. Since At = t, it follows that / = diag(^V aij, . . . , . a nj ), 
and / — ^4 = Li?. Then LB(x — v) = for every fixed point x of (6a). Because A is a consensus matrix, 1 is a 
simple eigenvalue of A, and At = 1. Then Ker(LS) = Ker(J — A) = Im(l). Notice that, by construction, rj = — Bv, 



L x, and Bt — 0. Consequently, (9 — rf) g Im(l), and statement (ii) follows. 



□ 



In Theorem 6.3 we assume that the measurements satisfy n g Im(B), or, equivalently, that the measurements are not 
affected by noise. While this assumption is justified by the fact that we only consider discretized measurements, in 
what follows we study the evolution of Algorithm 2 when the measurements are affected by noise. Let e(t) = r]—B9(t). 

Theorem 6.4 (Convergence of Algorithm 2 with noisy measurements) LetQ = (V,£) be a camera network, 
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and let Qd be its associated directed graph. Let r\ G 0™ be the vector of relative measurements, and let B be the 
incidence matrix ofQd- If A achieves average consensus over¥ p , then 

(i) there exists a finite time T < n such that the estimation error is constant, that is, e(t) = e(t + 1) for all t >T, 
and 

(ii) for all t>T, the estimation error satisfies 

e(t) = {^I-BY^A T L^ Vorth , 
where rjoith is the orthogonal projection of rj onto lu\{B)^- , and L is as in (7). 



PROOF: With the same notation as in the proof of Theorem 6.3, let T < n be the number of iterations needed 
for convergence of the network matrix A, that is, A T = n p ~ 2 tt T . Notice that BA t = for all t > T, so that 

t-1 T-1 

e(t) = v ~ A'as(O) - A T Lr/ = t] - ^ A T L V , 

for all t > T, and statement (i) follows. To show statement (ii), let r] — 77 par + ?7orth, where ?7p ar G lm(B), and 
??orth € Im^)^. From the linearity of (6a) and Theorem 6.3 we have 

T-1 

??par = B ^2 A T Ll] paI , 

T=0 

which concludes the proof. □ 



Theorem 6.4 characterize the performance of Algorithm 2 when the measurements are affected by noise. Notice 
that the estimation error can be minimized by properly choosing the network matrix A. We now conclude with a 
numerical example. 

Example 8 (Average computation and pose estimation via finite-field consensus) Consider a camera 
network with 4 cameras configured in a circle topology and network matrix 



A 



4 2 

110 4 

2 4 

12 3 



It can be verified that A achieves average consensus over F 5 in at most 3 iterations. In Fig. 5(a) we show that the 
network matrix A allows for the computation of the real-valued average of the agents initial states in 3 iterations 
(see Algorithm 1 and Section 4). 

Let Ak G jri° 24x1024 i, e the network matrix generated from A as A^ ~ A® A® A® A® A. The sparsity pattern of 
Ak is reported in Fig. 5(b). In Fig. 6 we validate our distributed pose estimation algorithm (see Section 6.2). □ 



7 Conclusion and Future Work 



In this paper we propose a distributed consensus algorithm for agents with limited memory, computation, and com- 
munication capabilities. Our approach is based on finite-fields, where agents states lie in a finite set, and operations 
are performed according to modular arithmetic. For our algorithm we identify necessary and sufficient convergence 
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(a) (b) 

Fig. 5. For the network described by the matrix A in Example 8, let 0, 1, 1, 1 be the agents initial states, respectively. Agents 
implement Algorithm 1. In Fig. 5(a) (above), we report the network disagreement as a function of time, where the disagreement 
at a given time equals the largest agents state minus the smallest agents state. Notice that the network achieves consensus in 
3 iterations (see Section 4). In Fig. 5(a) (below), we report the average computed by the first agent (solid black) as a function 
of time (see Section 6.1). Notice that the first agent, and hence every agent in the network, computes the average of the initial 
states (dashed red) at the third iteration. Fig. 5(b) shows the sparsity pattern of the matrix A k G p* 024 * 1024 in Example 8. 
The network defined by Ay. has 1024 nodes and 10 edges. 

conditions, and we characterize the convergence time. Additionally, we discuss several network design methods, 
and we propose some application scenarios. Our work proposes a novel class of consensus dynamics, which are 
advantageous in several applications, and it complements the existing literature on real-valued consensus. 

Through our analysis we show that finite-field consensus networks outperforms their real-valued counterpart in 
many aspects, including the convergence speed, robustness to communication errors, and agents requirements. These 
advantages come at the expenses of a more convoluted network design, which we identify as an interesting research 
direction. In particular, distributed network design algorithms, as well as gossip and asynchronous protocols should 
be investigated to broaden the applicability of finite-field consensus networks. Other theoretical research questions 
include characterizing the existence of finite-field consensus weights for a given interconnection structure and field 
characteristic, as well as the design of fastest finite-field consensus networks with fixed agents interconnection graph. 

APPENDIX 

Before proving Theorem 4.3, we recall the following fundamental results and facts in linear algebra. 

Theorem 7.1 (Primary decomposition theorem [25]) Let A : V — > V be a linear operator on some vector space 
V over some field ¥, and let p(s) — Yii=iPi( s ) be an annihilating polynomial for A with degree greater than 1, for 
some relatively prime polynomials pi, . . . ,p r . Then 

(i) Wi = Kei(pi(A)) is a A-invariant subspace for all i G {1, . . . , r}, 
(ii) V — Wi © W2 © • • • © W r , where © denote the direct sum operator, and 

(Hi) if Yii=i Pi( s ) — Pa(s) and ^-i *s the restriction of A to Wi, then pi is the characteristic polynomial of Ai. 

Recall that the order of a polynomial g G F[s], denoted by ord(g), is the smallest positive integer r such that g(s) 
divides s r — 1 over F, that is, the smallest positive integer r such that there exists q G F[s] satisfying s r — 1 = g(s)q(s). 



1G 





Time 

Fig. 6. For the camera network with 1024 cameras described by the matrix Ak in Example 8, this figure show the effectiveness 
of the our pose estimation algorithm in Algorithm 2. The figure above considers the noiseless case (Theorem 6.3). The dashed 
black line denotes the infinity norm of the estimation error, while solid red line corresponds to the (normalized) two norm of 
the difference between two consecutive estimation errors. We conclude that, in the absence of measurement noise, Algorithm 
2 converges in 3 iterations to a configuration compatible with the relative measurements. The figure below considers the 
case of noisy measurements (Theorem 6.4). Notice that, since the measurements are affected by noise, the estimation error 
remains nonzero (dashed black). However, the algorithm converges to a configuration with constant estimation error. Indeed, 
the difference between two consecutive estimation errors is zero after 3 iterations (solid red) . 

Theorem 7.2 (Order of a polynomial over finite field [24]) Let g £ ¥[s] be an irreducible polynomial satisfying 
g(0) 5^ and ord(g) = e. Consider f = g s , and t is the smallest integer such that p l > s, then ord(f) — ef>*. 

We are now ready to prove Theorem 4.3. 

PROOF of Theorem 4.3: Let Pa £ F p [s] be the characteristic polynomial of A, and notice that Pa can be 
written as 



P A {s) = det(sJ - A) = s h P(s), 
for some h £ N> , and P(s) £ F p [s] with P(0) ^ 0. 

(i) =>■ (ii) Due to Lemma 3.2, we have 1 £ o~ p (A). Thus, we factorize Pa in irreducible polynomials as 



(A-l) 



PA(s) = (s-l) k Y[Q J (s) m ^ 
where k, m,j £ N are given by the algebraic multiplicity of the corresponding eigenvalue. 

We start by showing that k = 1. Assume by contradiction that k > 1. Let W2 = Ker((I — A) k ), and let A2 be the 
restriction of A to W2 ■ Recall from [27] that the cycle structure of the transition graph Q 2 of A 2 is 



Cycles(£ 2 ) =d + Y,(P l ~ P l ^)Ci, 



(A-2) 



where the sums of cycles is just the corresponding union graph, and C\ denotes a unitary cycle, that is, a fixed point 
for A. From (A-2) it follows that, if k > 1, then the number of cycles in Q 2 is strictly greater than p. By Theorem 
4.1 we conclude that k = 1. 
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We now show that r = 0. Assume by contradiction that r > 0. Let W3 = Kcr(Qj(A) mj ), and let A3 be the restriction 
of A to W3. Then the cycle structure of the transition graph Q 3 of A 3 is 



^4 „deg(Q.,)i _ deg(Q.,)(i-l) 

Cycles(^) = Ci + £ ^ £ C /() 

i=i ^ 

where £j = ord(Q™ 3 ) > deg(Qj) > 1 from Theorem 7.2, and deg(-) denotes the degree of a polynomial. Since the 
graph structure of A is given by the product of the graphs associated with the irreducible factors of its characteristic 
polynomial [27], the number of cycles is greater than p whenever either k > 1 or r > (sec Example 4 and Fig. 2). 

(it) =>■ (i) Let Wi = Kct(A - I) = Im(l) and recall from Theorem 7.1 that Wi is A-invariant. Let V = [Vi 1] be 
an invertible matrix, where the columns of V\ are a basis for . Then we have 

A = v- x av = 

Since the eigenvalues of a matrix are not affected by similarity transformations, the characteristic polynomial of the 
matrix An is s n_1 , so that An is nilpotcnt. It follows that every vector in W x converges to the origin in at most 
n— 1 iterations, while vectors in Wi (consensus vectors) are fixed points for the matrix A. This concludes the proof. □ 
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